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CHAPTER I 
INTRODUCTION 


The heavy and diverse demands on aircraft for commercial, private, 
and military applications have made the task of the aircraft industry 
increasingly more complicated. Now millions of dollars must be in- 
vested over many years as teams of engineers carefully integrate the 
many systems involved in a wide-bodied jet or a fighter-bomber. Scale 
models are built and tested, and plans revised again and again as 
different possibilities are sorted through to find the few that are 
feasible and ultimately the one which best meets the demands for ap- 
plication and conditions for production. 

One of the systems which the design engineers must integrate 
into the total aircraft is the antenna system. A poorly functioning 
antenna on an automobile might result in the minor inconvenience of 
occasional loss of an AM radio signal. However, a similar loss of 
signal for an airplane performing an instrument landing could be fatal. 
Thus the number, type, and locations of antennas to be used is an 
important question, often made extremely difficult by the complicated 
nature of the solution of such electromagnetics problems. 

At low frequencies where the overall dimensions of the aircraft 
are on the order of several wavelengths or less, solutions to these 
problems can be obtained via the method of moments. And for extremely 
high frequencies, the classical methods of geometrical optics would 
provide a reasonable approximation. However, for a very wide range 
of frequencies, typically on the order of 10 y to Hertz both of 
these methods of solution break down. Moment method solutions become 
impractical as the number of unknown currents becomes too large to 
handle even on a large digital computer, and geometrical optics be- 
comes inaccurate as diffractions become significant. In the past 
the only practical method of predicting the performance of aircraft 
antennae in this frequency range was in the building of a scale model 
in order to experimentally determine its radiation patterns, a costly 
and time consuming process. 

However, recent advances in electromagnetic theory have made 
this task much simpler. Sommerfeld showed that the solution for a 
plane wave incident on a perfectly conducting wedge can be expressed 
as the sum of the geometrical optics wave (V*) and the diffracted 
wave (V B ) Til . And it was Pauli who made possible the expression of 
ML) as an asymptotic series [2] . Further work by many researchers 
such as Keller[3], and Hutchins and Kouyoumjian[4] lead to the Geo- 
metrical Theory of Diffraction (GTD), an asymptotic approximation 
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to the exact solution for diffractions from certain canonical forms, 
in the case where the dimensions of the scattering body are on the 
order of a wavelength or greater. 

Research in the GTD and its applications has been carried on 
at the ElectroScience Laboratory for almost two decades, resulting 
in increasingly more sophisticated solutions. Several years ago a 
computer code was developed, basically for far-field radiation pat- 
terns of antennae mounted on large commercial aircraft such as the 
Boeing 737, and delivered to NASA (Langley, Virginia) under Grant 
No. NGL 36-008-138. The aircraft design engineer was then provided 
with a relatively fast and efficient method of testing the performance 
of various antenna configurations. Since that time, improvements and 
modifications have been made to that code such that it can include 
greater application and increased efficiency. A major advancement 
has come about in the development of a new near field pattern solution 
for the Navy under Contract N00019-77-C-0299. 

The question to be dealt with here, dielectric slab scattering, 
was encountered when this near field code was applied to the study 
of small private aircraft [51. Although noteworthy achievements oc- 
curred in the modeling of various aircraft appendages such as fuel 
tanks and engines, this approach broke down in the modeling of a wind- 
shield, i.e., the only available solutions were for metallic bodies. 
For example. Figures 1 and 3 show two possible plate-cyl inder models 
that could be used to analyze the elevation plane pattern of a Cessna 
402B. The first assumes a transparent windshield (i.e., reflection 
coefficient of magnitude zero), and the second a totally reflective 
windshield (i.e., reflection coefficient of magnitude one). Compari- 
sons of the resulting computed patterns with experimental measure- 
ments provided by Melvin Gilreath at NASA (Langley, Virginia) are 
shown in Figures 2 and 4. Note that the experimental pattern has 
rather large oscillations in the region from the top backward 30° 
due to interactions between direct source fields and reflected fields 
from the windshield, which are not accurately predicted by either of 
the previous simulations. This results because the reflection coef- 
ficient of the windshield is not simply of magnitude zero or unity. 

In fact it is found to be generally a complex number dependent on 
the dielectric constant, thickness, frequency, and angle of incidence. 

The hypothesis presented here is that modifications may be made 
to geometrical optics and the GTD to provide accurate solutions for 
scattering from such dielectric layers. A derivation of the basic 
theory for reflection, transmission, and diffraction is given along 
with comparisons of predictions of this theory to known results and 
a summary. 
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(b) 


Figure 1. Computer simulation model for Cessna 402B elevation pattern 
with windshield assumed transparent and monopole antenna 
mounted forward of the cockpit. 

(a) top view, (b) side view. 
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Figure 2. Pattern resulting from Cessna 402B model with transparent 
windshield shown in Figure 1. 
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Figure 3. Computer simulation model for Cessna 402B elevation pattern 
with windshield assumed perfectly conducting and monopole 
antenna mounted forward of the cockpit. 

(a) top view, (b) side view. 
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Figure 4. Pattern resulting from Cessna 402B model with perfectly 
conducting windshield shown in Figure 3. 
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CHAPTER II 

TWO-DIMENSIONAL DIELECTRIC STRIP SCATTERING 


Consider the geometry illustrated in Figure 5 with a line source 
in the presence of a dielectric layer. The dielectric may be of any 
uniform material and thickness subject to the following restrictions: 
I) the diffractions may be assumed to emanate from a single point 
(Q); ?) energy leaving (Q) tangential to the layer due to a surface 
wave or transmission through the endface may be ignored; and 31 loss 
within the layer may be ignored. In light of our goal, a simulation 
model for an aircraft windshield, these restrictions are not overly 
severe. 


For a metallic half plane the total solution would be composed 
of incident, reflected, <sn«J diffracted fields; however, for the general 
discussion transmission must also be included. Thus the total electric 
field at an observation point (P) may be expressed as 


: tot = E V ♦ EV ♦ EV ♦ E d 


( 1 ) 


where u 1 , u r , and u^ are unit step functions shown here explicitly 
to emphasize the discontinuities in the geometrical optics (G.O.) 
solution. 


The individual terms of Equation (1) may be expressed, assuming 
a unit magnitude source, as 


E 1 (P) 


-jks. 

o 


Regions I and II 
Region III 


( 2 ) 


E r (P) 

E t (P) 



[ 0 


jr 

v. 'J i 


Region I 

Regions II and III 
Regions I and II 

Regions III 


(3) 

(4) 
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Regions I, II, and III 


. . - jks d 

E (P) = E 1 (Q)D Sjh ^ 

7 

(Here and throughout an e^ ot time dependence is assumed and supressed.) 
The "R", "T", and "D" are as yet undetermined coefficients of reflection, 
transmission, and diffraction, respectively. The task now is to de- 
termine expressions for these coefficients which when used in Equations 
(1) through (5) will yield the correct total field for scattering 
from a dielectric layer. In each case, the results will be stated 
for a perfectly conducting structure, then extended to include an 
arbitrary thin dielectric layer. 

A. Reflection 

If the half plane of Figure 5 is a perfectly conducting structure, 
the reflection coefficient will be of magnitude one with the sign 
depending on the polarization of the incident field. Thus, 

+1 magnetic line source 

R = "S (6) 

-1 electric line source 

for the case of perfect conductivity. 

For the more general case, consider a wave propagating in free 
space incident on a dielectric layer of thickness d and with a die- 
lectric constant e. as illustrated in Figure 6. At the front face 
the wave splits, a a portion of its energy being reflected and the rest 
being transmitted into the layer. At the opposite face the transmitted 
wave again splits, some energy being transmitted out of the dielectric 
and the rest reflecting back toward the front face. This process 
continues until all the original energy has left the dielectric, either 
as part of the total reflection or part of the total transmission. 

Thus the total reflection is an infinite sum of waves each of 
which can be related back to the incident wave provided that a plane 
wave propagation approximation may be made. Letting R, represent 
the reflection coefficient for the initial external reflection, R ? 
represent subsequent internal reflection coefficients, T, and T ? repre- 
sent transmission coefficients into and out of the layer, respectively, 
represent the phase delay associated with the field in a single 
crossing of the slab, and P be a term to account for the difference 
in path length to the observer for different rays leaving the layer 
in the case of oblique incidence as illustrated in Figure 7, the total 
reflected field may be expressed as 
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E tot = {R 1 + T l T 2 l ( R 2 ) 2n " 3 ( p d )2n " 2 ( p a )n_1 } El • (7) 

n=2 

Use of P. as in the above equation assumes plane wave propagation 
within the dielectric. And this usage of P assumes that the rays 
between the slab and the observation point are parallel. These as- 
sumptions are satisfied provided that neither the source nor receiver 
is too close to the slab. These assumptions and some basic trigonom- 
etry yield 


P d =e- 3k ' 8, (8) 

J k2Jlsin0^si n0 . 

P =e 11 0) 

a 


where £,=d/cose t and (k) and (k 1 ) are the propagation constants for 
free space arid the dielectric medium, respectively. Note that k'=k/e^, 
where e r = e d /e 0 . 


The coefficients of single reflection and transmission are found 
via continuity of the tangential components of the E and H fields, 
and are given by[6] 


R KM) 


e^cos© .- 

I 1 nJ 

. 2 

e r -sin 0 

e„cose .+ 

1 1 S 

e -sin 2 © . 

t- r 


(10a) 


(l/ejcos0.- / (l/ej-sin 2 e t 

R 2(M) ' — <H«) 
(l/e r )cos 0 t + J(l/e r )-sin 0 t 


for the case of a magnetic line source (E field parallel to the plane 
of incidence), and 


‘1(E) 


cose i - w 


e r -sin 2 0 i 


2(E) 


cos0 i + i er- sin2a i 

cos0 t - J(l/e r )-sin 2 0 t 
cos0 t + n /(1\) -sin 2 0 t 


(10b) 


(lib) 


•for the case of an electric line source (F field perpendicular to 
the plane of incidence). In both cases one finds 
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( 12 ) 


T 2 = R 2 + 1 

Op 0 "j 

sino t = sine i 4r^;. 


(13) 

(14) 

(15) 


Through the evaluation of the geometrical series and noting that 
R 2 --Rp Equation (7) may be simplified to the desired total reflection 
coefficient: 


R (M,E) 


R l(M t E)^- p d p a) 
1_R 1(M,E) P d P a 


(16) 


B. Transmission 


Referring again to Figure 5, one finds for a perfectly conduct- 
ing half plane that T = 0. For the dielectric case, however, the 
transmitted field is once again an infinite sum and may be expressed 
as 


E 


tot 


■ (Wt l =1 (p d ) 2 n ~V/'V 2 ) 2n " 2 


JE 1 


(17) 


where E 1 is the field that would have been incident on the receiver 
in the absence of the dielectric layer. Here P, is a term to cause 
the transmission coefficient to act as a simple modification of the 
incident field as shown in Figure 8(a). Note that without this term 
the phase is referred to point (B) in Figure 8(b). This term is given 
by 

Jk*cos(9 r e t ) . ( 18 ) 

H t e 


As before in the case of reflection, much simplification is pos- 
sible. After some manipulation one arrives at the desired total trans- 
mission coefficient: 


T (M,E) 


( 1 ~ R UM, E)) P d P t 

1-R l(M,E) P d P a 


(19) 
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C. Diffraction 


Consider the geometry illustrated in Figure 9 which shows a wave 
incident on a perfectly conducting wedge, one of the canonical forms 
for which the Geometrical Theory of Diffraction (GTD) provides an 
asymptotic solution. The diffraction coefficient for the 2D case 
is generally given by 


D = 


r 



magnetic line source 
electric line source, 


( 20 ) 


where the "h" and "s" subscripts refer to acoustically "hard" (E field 
perpendicular to diffracting edge) and "soft" (E field parallel to 
diffracting edge) diffractions, respectively. Reference [7,81 gives 
these coefficients for a half plane (n=2 in Figure 9) as 

_j I f ■) 

D £-== 1 - — /- fer P- 1 * 1 . (2D 

■’ 2/2nk sin3 Q cos (^r~j cos j 


INCIDENT REFLECTION 
WAVE SHADOW BOUNDARY 



Figure 9. Geometry for perfectly conducting wedge diffraction. 
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where 


F(X) = 2jJX e J ' X / 

77 



ch . 


Often the notation 

D s 3 h ^ ( l >s ^ * ^ = + D(tj)+<Jj 1 ) 


( 22 ) 


where . jt 

= zg — 1 £lkU(*ii^.L 

277k sine 0 cos ^£|l) 

is used. This form shows explicitly how the GTD coefficients com- 
pensate for each of the two discontinuities in the G.O. solution. 

The D ( (})-cp ’ ) term is associated with the incident shadow boundary and 
the Df<j>+<|>') term is associated with the reflection shadow boundary. 

Note that in the application of Equation (22) one would use the 
(-) sign for the case of an electric line source and the (+) sign 
for the case of a magnetic line source. Thus the "R" of Equation 
(6) may be incorporated here, and the diffracted field written as 

_J ’ ks d 

E d (P) = E 1 (Q) [D((j)-()) 1 )+R D(<f>+<t>')J — z. — . (23) 

^d 

Note that for the perfectly conducting case the discontinuity at the 
incident shadow boundary is unity times the incident field near the 
boundary and the coefficient of D( c|>-(J) ' ) is unity; likewise, the dis- 
continuity at the reflection shadow boundary is R times the reflected 
field near that boundary and the coefficient of D ( cj>+d> 1 ) is R. 

As the subject of attention is now shifted to diffractions from 
dielectrics, the discussion must be limited to flat layers rather 
than wedges since uniform thickness was a requirement for derivation 
of the reflection and transmission coefficients. As noted above, 
there is a correlation between the magnitude of the discontinuity 
at each shadow boundary and the coefficient of the corresponding Dkh+iji 1 ) 
term. In the dielectric case the discontinuity at the reflection 
shadow boundary is again related to R, although R is no longer simply 
-1; at the incident shadow boundary the discontinuity is O-T) times 
the incident field, which reduces to unity for the metallic case. 
Therefore, it is conjectured that the diffraction coefficient for 
a finite dielectric layer is given by 

D = r(l-T)D( ( j)-(l) l )+R D ( (f>+<)) 1 ) ] . (24) 
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Ilia values of "R" and "T" are given by Equations (16) and (19) with 
0. being the angle between the ray incident on the diffraction edge 
and the normal to the reflecting face. 

Equations (16), (19), and (24), when used in Equations (1) through 
(5) give the desired solution for the problem shown in Figure 5. This 
has been demonstrated experimentally as will be discussed in Chapter 
IV. However, before presenting this confirmation, the theory will 
be expanded to the more general three dimensional case. 
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CHAPTER III 

THREE DIMENSIONAL DIELECTRIC SLAB SCATTERING 


In the previous chapter coefficients were developed for the var- 
ious terms in the G.O.-GTD solution for two dimensional scattering 
from a finite dielectric layer. These coefficients are still valid 
in the more general case, however, care must be taken in their ap- 
plication. In treating 3D problems one must keep in mind two facts: 

1) spherical waves are being considered rather than cylindrical waves, 
and ?.) the field can be arbitrarily polarized rather than merely that 
of a simple electric or magnetic line source. The first of these 
facts means that the spreading factors become more complicated than 
e~ JKS /-/s; the second means one must find a suitable coordinate system 
in which to express the resulting fields, transforming coordinates 
as necessary to apply the coefficients of the previous chapter. 


It will be convenient at^this point to introduce the "ray fixed 
coordinate system." Letting n be the unit vector normal to the slab, 
I be,,the incident unit vector from the source to a point on the slab, 
and R be the reflection unit vector from the point of reflection to 
the observation point, one may define the following: 

A A A 

u A = n x I, 

A j A /\ 

U„ = I X u A , 

/\ V* /\ /S 

U„ = R X U A , 


where (•*•) and (") indicate vectors perpendicular and parallel, re- 
jjpectiyely, to the ordinary plane of incidence, the plane containing 
I and n. These unit vectors form the orthogonal bases for the ray 
fixed coordinate system as illustrated in_ N Figure 10. The field of 
a ray traveling in the direction of I or R is. now completely speci- 
fied by its components E^and E lt . Note that E^ corresponds to the 
2D case of an electric line source and E, t corresponds to the 2D case 
of a magnetic line source. In the ray fixed system one may utilize 
the coefficients of the previous chapter and express the reflected 
and transmitted fields in the following manner: 




R (M) 

0 


0 

R (E)- 



f r (s) 


(28) 
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0 


(29) 



T 

0 


(M) 


E) 



f t (s) 


The angle of incidence to be used in computing the coefficients is 
given by 


cos 


( -I *n) 


(30) 


The spreading factor for the transmitted field is the same as 
that for an incident field, 


f 


t 



(31) 


For the reflected field the 
plicated, given by T8 1 


r r 


f r (s) 


p l p 2 




spreading' factor is slightly more com- 


e 


-jks 


(32) 


The ray fixed coordinate system, although ideal for reflection 
and transmission problems, is not suitable for the diffraction solu- 
tion. The diffraction coefficients presented in the previous chapter 
are valid only when the incident field is expressed by components 
parallel and perpendicular to the edge fixed plane of incidence, the 
plane containing the incident ray and the diffracting edge; and the 
diffracted field expressed by components parallel and perpendicular 
to the edge fixed plane of diffraction, which contains the diffracted 
ray and the edge. Thus the "edge fixed coordinate system" is intro- 
duced and illustrated in Figure^ll. Letting e be the unit vector par- 
allel to the diffracting edge, I be the incident unit vector from 
the source to the diffraction point (CL), and 0 be the diffraction 
unit vector from (Q E ) to the observation point (P), one may define 
the following: 


$' = -e x I, 

(33) 

Pi = $’ x I, 

(34) 

$ = e x D, and 

/V A A 

(35) 

3 = <h x D. 

1 36) 


These vectors form the two orthonormal bases of the edge fixed co- 
ordinate system. Note that the ordinary plane of incidence inter- 
sects the edge fixed plane of incidence along the incident ray, and 
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the edge fixed plane of diffraction along the diffracted ray. The 
relationship between these planes, and thus the relationship between 
the ray fixed system and the edg i fixed system is depicted in Figure 
12. It can be shown [81 that the angles a and a., defined as shown 
in Figure 12 are related by «. = -a. 

It is apparent that the incident field may be expressed in the 
edge fixed coordinate system as 

Eg, = E, 1 , cosa - e] sina (37a) 

' o 

eJ, = E, 1 , sina + e] cosa . (37b) 

These expressions may be written more compactly in matrix notation 
a j 


E 


i 

e.f . 


Tf-ajE 1 


where 


(38) 


T / „ cosa -sina' 

T( -a; - 

.sina cosaj 

Likewise the reflected field is transformed between systems 


(39) 

by 


E^ = E^ cosa + Ex sina (40a) 

p o 

E^ = -En sina + E^ cosa . (40b) 

This may be expressed as 


E e.f. - T <“> Er • 

where 


(41) 


T(a) 


cosa sina 
-sina cosa 


(42) 


In the 
in a manner 


edge fixed system, the 3D diffracted field may be expressed 
analogous to the 2D case by 
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Figure 12. Comparison of the ray fixed coordinate system and the 

edge fixed coordinate system. The incident and reflected 
rays are perpendicular to the page and directed outward. 




whf»r<> the matrices A and B serve the same function of scaling the 
'Kijti'b'l terms to the proper complex value of the discontinuity as 
did (R) arid H-Tl in the 2D case. A is the. polarization matrix as- 
sociated with the discontinuity in the G.O. field at the incident 
shadow boundary and R the polarization matrix associated with the 
discontinuity at the reflection shadow boundary. However, These dis- 
continuities must he expressed in terms of the unit vectors § , $. 

On the lit side of the incident shadow boundary one finds (see 
Figure 1 1 ) 



On the side not directly illuminated by the source, the transmitted 
field is present. To express this one must start with E in the edge 
fixed system, transform to the ray fixed system in which the trans- 
mission coefficients are valid, perform the transmission operation, 
then transform back to the edge fixed system. Thus the field in this 
region is given by 



So the total discontinuity is given by 



( 45 ) 



such that 

2 2 

-1 + T^cos a + T^sin a (T^ - T^)sina cosa 

2 ? 

(T(m) " T(£))sin<x cosa -1 + T^sin a + T^cos a 

( 46 ) 
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Likewise at the reflection shadow boundary one notes that 



Therefore the discontinuity is given by 



such that 

2 2 

R^cos a- R^sin a (R^ + R^)sina cosa 
B = ' 2 ? (48) 

- (R (M) +R ( E ) )sina C0Sa - R (M) sin a + R (E) C0S ' a J ’ 


For computational purposes, the 3D diffracted field is written 
as 



where 

D a = (l“ T (|vi) cos2a " T (E) s ' ir|2a )D((!>-<))' )-(R( M )C° s2 a-R( E ^sin 2 a)D((!)+())' ) 

(50) 

D b = ( ~ ( ^" ( M ) ~T ( e j) s i nacosa) D ( <|)-<|> 1 ) - ( (R +R ^ ^ j)s i nacosa) D ( d>+<j> ' ) 

(51) 


° c = (~(^(m)~T( E ) )si nctcosa) D ( 4>-4> ’ ) + ( (R^^+R^ )sinacosa)D((j)+d) 1 ) 

(52) 

D d = ^ 1-T (M) si,1 ^ ot ~ T (E) COS ^ a ^^~ < ^ 1 ) + ( R (|v]) sir|2a '' R (E) cos ^ a ) D ( < l )+ ^' ) • 

(53) 
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The total 3D solution for scattering from a dielectric slab is 
now complete, given by the sum of Equations (28), (29), and (49). 

In the next chapter, computed patterns are compared to experimental 
results for confirmation of the theory. 
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CHAPTER IV 
RESULTS 


A. Coefficients 


»j.p»05)UCIB1UTY OF THE 
original page IS POOR 


Plots showing the variation of the reflection and transmission 
coefficients of Chapter II with respect to the dielectric layer thick- 
ness are shown in Figures 13(a) and 1. 3 f b ) . respectively. These re- 
sults are for a typical dielectric (polystyrene) illuminated by a 
plane wave with normal incidence (Q.=Cr in Figure 6) In the eval- 
uation of the reflection coefficient, one would expect the first two 
terms in the series to dominate. These two terms will differ in phase 
by the propagation length of two transversals of the layer plus a 
180° phase shift due to the R ? reflection off the back face Thus 
a maximum is anticipated at tne points where D/X = (2n+l)/4 with 
n=l , 2 , 3 , • * * , whicli is the case of the propagation phase delay p 1 y s 
the reflection phase change equaling an integral multiple of 360 , 
i.e., the first two terms of the series adding in phase. A minimum 
would be expected when these terms add out of phase, which occurs 
at D/X = n/2. Both predictions match the results shown in Figure 
13. 


In the transmission case, the second term has two R 0 reflections, 
which add a 360° phase change. Thus the phase difference between 
the two dominant terms is effectively caused by the propagation of 
two layer transversals. Thus maxima occur at D/X = n/2 and minima 
at D/X = (2n+U/4. These again match the computed results shown in 
Figure 13. 

B . Dielectric Layer Illuminated by a Rectangular Horn 

For purposes of an initial testing of the theory in a real world 
situation it is applied to the problem shown in Figure 14, where the 
pattern is measured in the E plane of the horn with the phase origin 
at the center of the horn aperture. A corresponding simulation model 
made use of the above theory for scattering from the polyethylene. 

The horn was treated as a set of magnetic line sources (one at the 
apex and one at each of the diffracting edges) as in the corner re- 
flector model presented in Reference [91. The first measurements were 
taken with the horn alone, and the results at various frequencies 
are shown in Figures 15 through 19. Since only the main beam of the 
horn will interact with the dielectric, and only these interactions 
are of interest, extreme care was not taken in the setup of the horn. 
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Figure 13a. Plot of reflection coefficient vs. slab thickness for normal incidence on a 
polystyrene slab (e r =2.55) with the E field normal to the plane of incidence. 
Note that "D" is measured in wavelengths within the dielectric. 
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Figure 13b. Plot of transmission coefficient vs. slab thickness for normal incidence on a 
polystyrene slab (e =2.55) with the E field normal to the plane of incidence. 

D is measured in wavelengths within the dielectric. 
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Figure 14. Polyethylene geometry used to test scattering theory. 

Measurements are taken in the E plane of the horn. 
Polyethylene thickness is .1285" and e =2.25 
at 10 GHz. A standard gain horn manufactured 
by Narda Corp. is used. 
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This accounts for much of the difficulties in the back region of the 
horn, especially 0 to 60 and 300 to 360 . Most notably, the meas- 
ured patterns are not symmetrical. However, in all cases the main 
beam is in good agreement with the computed results. This, as 
will be seen, is sufficient for the testing of the theory. 

Results for the entire geometry illustrated in Figure 14 are 
shown in Figures 20-24 for various frequencies. Again difficulties 
are observed in the region behind the horn, 0° to 60° and 300 u to 
360°; however, these can be attributed to the experimental setup rather 
than the dielectric solution, which is the subject of investigation 
at present. For confirmation of the theory, one should look to the 
two areas of the pattern most strongly affected by the dielectric, 
i.e., the area of direct transmission around <t>=180 u and the area of 
direct reflection around <j)=90°. In these two areas one finds the 
results quite acceptable. The shape of the main beam, which is de- 
termined totally by transmitted and diffracted fields, is at all fre- 
quencies found to be in good agreement. In addition, the reflection 
region shows good results. Here ripple is produced by interactions 
between the horn sidelobes and scattering from the dielectric. The 
amplitude and frequency of this ripple agree reasonably well with 
the measured ripple in each case. 

Figure 25 shows a second geometry used for verification of the 
dielectric scattering theory. Here the slab is much thicker and of 
slightly higher dielectric constant. Also it is closer to the horn 
which will reduce the spacial frequency of pattern ripple, making 
it easier to compare results. Calculated and measured patterns shown 
in Figures 26-28 for various frequencies agree well in overall shape, 
especially in the main beam and reflection regions. 

Figure 29 shows another geometry with the same slab as the pre- 
vious problem placed perpendicular to the horn axis. In th i g case 
the region of primary interest is the main beam around (|>=180 , deter- 
mined completely by transmitted and diffracted fields. The calculated 
and measured patterns illustrated in Figures 30-34 show reasonably 
good agreement in each case. 

C. Dipole Illuminating a Dielectric Covered Ground Plane 

The geometry illustrated in Figure 35 shows a test of the theory 
using a dipole source rather than a horn, and with the dielectric 
slab backed by a ground plane. Here the calculated results are found 
by a re-derivation of the coefficients of Chapter II setting the coef- 
ficient for reflections from the second face equal to -1, and the 
experimental measurements were obtained by measuring the pattern in 
the planes defined by the z axis and each of the dashed lines of Figure 
35. Results are shown in Figures 36 and 37 for the principal planes 
and Figures 38 and 39 for the two polarizations of the diagonal pattern. 
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Figure 20. Pattern for scattering from polyethylene strip of 
thickness .1285" at 8.0 GHz. 
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Figure 21. Pattern for scattering from polyethylene strip at 9.0 GHz. 
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Figure 22. Pattern for scattering from polyethylene strip at 10.0 GHz. 
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Figure 23. Pattern for scattering from polyethylene strip at 11.0 GHz. 
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Figure 24. Pattern for scattering from polyethylene strip at 11.9 GHz 





Figure 25. First polystyrene geometry. Slab thickness is .2715" 
and e =2.55. Measurements taken in the E plane of the 
horn. Results shown in Figures 26-28. A standard 
gain horn manufactured by Narda Corp. is used. 
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Figure 26. Pattern for polystyrene geometry shown in Figure 25 at 8. GHz. 
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Figure 28. Pattern for polystyrene geometry shown in Figure 25 at 11.5 GHz 




Figure 29. Second polystyrene geometry used especially for 
investigation of transmission. Again slab 
thickness is .2715" and e =2.55 with measurements 
taken in the E plane of tne horn. Results shown 
in Figures 30-34. A standard gain horn manufactured 
by Narda Corp. is used. 
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Figure 30. Pattern for polystyrene geometry shown in Figure 29 at 8. GHz. 



Figure 31. Pattern for polystyrene geometry shown in Figure 29 at 9. GHz. 
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Figure 33. Pattern for polyethylene geometry shown in Figure 29 at 11. GHz. 
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Figure 34. Pattern for polyethylene geometry shown in Figure ?9 at 11 .5 GHz. 
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Figure 35. A polystyrene slab which is .2715" thick and e =2.55 is 
mounted flush on a square ground plane. The structure 
is illuminated by a dipole mounted 5.15" above the cen- 
ter of the plate. Using normal spherical coordinates, 
the patterns are measured by varying 6 with <j> a constant 
Note that the z axis is into the page. 
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Figure 36. E pi 



pattern for geometry illustrated in Figure 35. 
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Figure 37. H plane pattern for geometry illustrated in Figure 35. 









In order to illustrate the effect of the dielectric layer, a 
comparison is made between these patterns and corresponding patterns 
obtained from the same geometry in which the dielectric layer is re- 
moved in Figures 40-43. It is apparent that the dielectric does sig- 
nificantly affect the pattern making this geometry a valid test of 
the dielectric theory. 

D. Comparison with Moment Method 

Further confirmation is obtained, as well as a study of appli- 
cability by comparison of results of this theory with Moment Method 
results [10,111 using the geometry shown in Figure 44. The moment 
method results, as presented here, are assumed to be the exact solu- 
tions to these problems. Working with first an electric and then 
a magnetic line source, the effects of varying source to strip distance 
and angle of incidence are studied. Results for varying the source 
range are shown in Figures 45 and 46. Here it is noted that our re- 
sults are accurate at ranges of greater than a wavelength, but begin 
to break down inside that range. Results for varying the angle of 
incidence are shown in Figures 47 and 48. It is apparent here that 
the solutions are accurate to at least 40° off normal incidence and 
can perhaps be extended to 60° off normal incidence. 
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Figure 41. Comparison of H plane patterns of qeometry illustrated in Figure 35 
with (solid line) and without (dashed line) the polystyrene slab. 
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Figure 43. Comparison of <f> components of diagonal patterns of geometry illustrated in Figure 
35 with (solid line) and without (dashed line) the polystyrene slab. 

(Receiving horn in E plane . ) 
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Figure 44. Geometry used for GTD-Moment Method comparison. 

Note that 4> and p are varied for an electric 
and a magnetic line source. 
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(a) ^ 


Figure 45. GTD-MM comparison with electric line source and varying p using geometry shown 
in Figure 44 with L=2.X, D=.05X, <b s = 0. (a) P S = 2.X. (b) p^-l.X. 



(c) fd) 

Figure *5. GTD-MM comparison with electric line source and varying p , using geometry shown 
in Figure 44 with L=?.X , D=.05\, $ =0. 1 c) p =.5X. 1 d~ 1 p =.2^X. 
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Figure 46. GTD-MM comparison with magnetic line source and varying p s , using geometry shown 
in Figure 44 with L=1.A, D=.I\ , ($l= 0. ^ a ) P s = 2-A. (b) p s =l-A. 
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Figure A F. 


GTO-MM comparison with magnetic 
in Figure 44 with L=l.A, D=.U, 


line source and varying p usinq geometry shov/n 
4> s =0. (c) P S =-5X. (d) p*=.25X.~ 





(a) 


'10 


Figure 47. GTD-MM comparison with electric line source and carrying <j> , using geometry shov/n 
in Figure 44 with L=2.A, D=.05A, p s =2.A. fa) 4> s =20 . . (bi 4> S =40 Q . 
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Fiqurc 48. GTD-MM comparison with magnetic line source and varying <j> , using geometry shown 
in Figure 44 with L=1.X, D=.lX 9 (d= 2.X. (a) d> s = 20. . (bj ^-40. . 




Figure 48. GTD-MM comparison with magnetic line source and varying $ , usina geometry shown 
in Figure 44 with L=1 , D=. 1>, , p $ =2. . (c) <j> s =60 . ( d ) s ^ s ~80 . . 



CHAPTER V 
CONCLUSIONS 


This thesis has demonstrated that geometrical optics and the 
Geometrical Theory of Diffraction can be modified to include scat- 
tering from thin dielectric slabs. Thus it will soon be possible 
to model structures such as private aircraft in which non-metal lie 
parts play a significant role. This approach will provide efficient 
computation of a volumetric antenna pattern and is expected to give 
reasonable accuracy. 

The limitations on the application of this theory are not severe, 
but must be noted. As shown in the Moment Method comparisons, the 
source must be kept at least a wavelength or so away from the dielectric 
to keep within the plane wave approximation of Chapter II. The angle 
of incidence, i.e., the angle between the incident ray and the normal 
to the dielectric layer, must be kept to not much greater than 60 ° 
or so to avoid exciting a surface wave, which was assumed not signif- 
icant in Chapter II. An additional limitation is on the thickness 
of the dielectric. The theory was developed on the assumption of 
the thickness being less than one tenth wavelength; however, accurate 
results have been obtained for dielectrics with thickness of nearly 
one half wavelength in free space based on comparisons with measure- 
ments. 

The question of diffractions from the junction of a dielectric 
slab and a metal plate or two dielectric slabs remains unanswered, 
but it is expected that a solution is available by the same methods 
presented here. Derivation of this would then provide the total theory 
needed for full incorporation in an aircraft antenna code. 
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